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We study the linear response of a coherently driven polariton fluid in the pump-only configuration 
scattering against a point-like defect and evaluate analytically the drag force exerted by the fluid on 
the defect. When the system is excited near the bottom of the lower polariton dispersion, the sign 
of the interaction-renormalised pump detuning classifies the collective excitation spectra in three 
different categories [C. Ciuti and I. Carusotto, physica status solidi (b) 242, 2224 (2005)]: linear for 
zero, diffusive-like for positive, and gapped for negative detuning. We show that both cases of zero 
and positive detuning share a qualitatively similar crossover of the drag force from the subsonic to 
the supersonic regime as a function of the fluid velocity, with a critical velocity given by the speed of 
sound found for the linear regime. In contrast, for gapped spectra, we find that the critical velocity 
exceeds the speed of sound. In all cases, the residual drag force in the subcritical regime depends 
on the polariton lifetime only. Also, well below the critical velocity, the drag force varies linearly 
with the polariton lifetime, in agreement with previous work [E. Cancellieri et al, Phys. Rev. B 
82, 224512 (2010)], where the drag was determined numerically for a finite-size defect. 
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I. INTRODUCTION 

Out of equilibrium quantum fluids such as polaritons 
in semiconductor microcavities are being the subject of 
an intensive study. Microcavity polaritons, the quasipar- 
ticles resulting from the strong coupling of cavity photons 
and quantum well excitons, have the prerogative of be- 
ing easy to both manipulate, via an external laser, and 
detect, via the light escaping from the cavity pQ. In par- 
ticular, resonant excitation allows the accurate tuning 
of the fluid properties, such as its density and current. 
However, the polariton lifetime being finite establishes 
the system as intrinsically out of equilibrium: An exter- 
nal pump is needed to continuously replenish the cavity 
of polaritons, that quickly, on a scale of tens of picosec- 
onds, escape. 

Recently, the superfluid properties of a resonantly 
pumped polariton quantum fluid in the pump-only con- 
figuration — i.e., where no other states aside the pump 
one are occupied by, e.g., parametric scattering — have 
been actively investigated both experimentally and the- 
oretically 2-9 . This pumping scheme, differently from 
other cases, such as the resonant optical parametric oscil- 
lator regime and the non-resonant pumping scheme, cre- 
ates a polariton fluid that, inside the pump spot, is not 
characterised by a free phase. On the contrary, the phase 
of the pump state is locked to the one of the external 
pumping laser. Nevertheless, it has been predicted [3] 
and observed [4] that scattering can be suppressed below 
a critical velocity, where the system displays superfluid 
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behaviour, similarly to what has been predicted by the 
Landau criterion for equilibrium superfluid condensates. 
Further, a fixed phase clearly prevents the formation of 
phase dislocations, such as vortices and solitons. For 
this reason, it has been suggested [S] and experimentally 
realised [7j that the defect can be located just outside 
the pump spot, where the hydrodynamic nucleation of 
vortices, vortex-antivortex pairs, arrays of vortices, and 
solitons can be observed when the fluid collides with the 
extended defect. Similarly, nucleation of vortices in the 
wake of the obstacle has been observed in pulsed experi- 
ments 0GI]. 

In a conservative quantum liquid flowing past a small 
defect, the Landau criterion for superfluidity links the on- 
set of dissipation at a critical fluid velocity with the shape 
of the fluid collective excitation spectrum [10]. In par- 
ticular, for weakly interacting Bose gases, the dispersion 
of the low-energy excitation modes being linear implies 
that the critical velocity for superflow coincides with the 
speed of sound c s . Clearly, this is strictly correct only 
for vanishingly small perturbations |11) . while for a de- 
fect with finite size and strength, the critical velocity can 
be smaller than c s [HI [T3] . 

However, even for perturbatively weak defects, in out- 
of-equilibrium systems, where the spectrum of excita- 
tions is complex, the validity of the Landau criterion has 
to be questioned HH [15]. In the particular case of 
coherently driven polaritons in the pump-only configura- 
tion, it has been predicted [5J |3] , and later observed [1] , 
that scattering is suppressed at either strong enough 
pump powers or small enough flow velocities. Yet, on 
a closer scrutiny, it has been shown that, despite the 
apparent validity of the Landau criterion, the system al- 
ways experiences a residual drag force even in the limit of 
asymptotically large densities [5] or small velocities. This 
result has been proven by numerically solving the Gross- 
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Pitaevskii equation describing the resonantly-driven po- 
lariton system in presence of a non-perturbative extended 
defect. Here, the drag force exerted by the defect on the 
fluid has been shown to display a smooth crossover from 
the subsonic to the supersonic regime, similarly to what 
it has been found in the case of non-resonantly pumped 
polaritons [15 . In this work, we find an even richer phe- 
nomenology for the dependence of the drag force on the 
fluid velocity and two different kinds of crossovers from 
the sub- to the supercritical regime. Further, we show 
that the origin of the residual drag force, which, in agree- 
ment with Ref. }5] , lies in the polariton lifetime only, can 
be demonstrated even within a linear response approxi- 
mation. 

More specifically, in this work, we apply the linear 
response theory to analytically evaluate the drag force 
exerted by the coherently driven polariton fluid in the 
pump-only configuration on a point-like defect. To sim- 
plify the formalism, we restrict our analysis to the case 
of resonant pumping close to the bottom of the lower 
polariton dispersion, where the dispersion is quadratic. 
Here, the properties of the collective excitation spectrum 
have been shown to be uniquely determined by three pa- 
rameters only [3j: the fluid velocity v p , the interaction- 
renormalised pump detuning A p , and the polariton life- 
time k. In particular, the sign of the detuning A p deter- 
mines three qualitatively different types of spectra: lin- 
ear for A p — 0, diffusive-like for A p > 0, and gapped for 
Ap < 0. 

For both cases of linear and diffusive spectra, we find 
a qualitatively similar behaviour of the drag force as a 
function of the fluid velocity v p : In particular, the drag 
displays a crossover from a subsonic or supcrfluid regime 
- characterised by the absence of quasiparticle excita- 
tions — to a supersonic regime — where Chcrcnkov-like 
waves are generated by the defect and propagate into the 
fluid. The crossover becomes sharper for increasing po- 
lariton lifetimes 1/k and displays the typical threshold 
behaviour for k — ¥ with a critical velocity given by the 
speed of sound of linear regime, v c — c s , exactly as for 
weakly interacting equilibrium superfluids (in the case 
of perturbatively weak defects). This behaviour is simi- 
lar to the one predicted for polariton superfluids excited 
non-resonantly |15j . where the spectrum in that case is 
diffusive- like. 

However, for gapped spectra at A p < 0, we find that 
the critical velocity governing the drag crossover exceeds 
the speed of sound, v c > c s , and we determine an an- 
alytical expression of v c as a function of the detuning 
Ap. Further, for n — > 0, the drag has a threshold-like 
behaviour qualitatively different from the one of weakly 
interacting equilibrium superfluids, with the drag jump- 
ing discontinuously from zero to a finite value at v p = v c . 

We evaluate the drag as a function of the polariton 
lifetime k and find for all three cases that: In the super- 
critical regime, v p > v c , the lifetime tends to suppress 
the propagation of the Cherenkov waves away from the 
defect and therefore to suppress the drag. Instead, well in 



the subcritical regime, v p <C v c , we find that the residual 
drag goes linearly to zero with the polariton lifetime n, in 
agreement to what it was found in Ref. [5] , by making use 
of a non-perturbative numerical analysis for a finite size 
defect. Similarly to Ref. [5], here, we do also find that the 
residual drag in the subcritical regime can be explained 
in terms of an asymmetric perturbation induced in the 
fluid by the defect in the direction of the fluid velocity. 

This paper is structured as follows: In Sec.|H]we briefly 
introduce the linear response approximation. We classify 
the three types of collective excitation spectra in the sim- 
plified case of excitation close to the bottom of the lower 



polariton dispersion in Sec. II A In Sec. Ill we derive 
the drag force and characterise the crossover from the 
subsonic to the supersonic regime in the three cases of 
zero, positive and negative detuning. In this section, we 
also evaluate the drag as a function of the polariton life- 
time, interpreting therefore the results of Ref. [5j- Brief 
conclusions are drawn is Sec. IIVI 



II. LINEAR RESPONSE 

The description of cavity polaritons resonantly excited 
by an external laser is usually formulated in terms of 
a classical non-linear Schrodinger equation (or Gross- 
Pitaevskii equation) [T^] for the lower polariton (LP) field 
V>lpM) (H=l): 



idtipLP = [w L p(-iV) - in + V(r) + g\^ LP \ 2 ]ilj 



LP 



+ F(r,t). (1) 

The LP dispersion is expressed in terms of the photon 



2m c 



and exciton lo x energies, the photon 



mass mc, and the Rabi splitting pQ: 



LO LP (k) = - [ WC (k)+4] 



[w c (k) - ' 



'XI 



(2) 



Because polaritons continuously decay at a rate k, the 
cavity is replenished by a continuous wave resonant pump 
F(r, t) at a wavevector k p (we will later assume k p di- 
rected along the ^-direction, k p = (k p ,0)) and frequency 



T(r,t) = f p e 



-UJpt) 



(3) 



Note that, as discussed in appendix |A| Eq. ([TJ is a sim- 
plified description of the polariton system: This implies 
that the interaction non-linearities are small enough not 
to mix the lower and upper polariton branches. More- 
over, starting from a formulation in terms of coupled 
exciton and photon fields, the polariton lifetime would 
be momentum dependent and, similarly, the polariton- 
polariton interaction strength g is not contact-like as in- 
stead assumed in Eq. ([!]). However, as shown in ap- 
pendix^ these simplifications, do not affect our results 
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qualitatively, rather, allow to write them in terms of sim- 
pler expressions. Further, we have checked that, when- 
ever the system is excited near the bottom of the lower 
polariton dispersion, the results for the drag force re- 
ported in Sec. Ill coincide with the ones obtained by 



using an exact photon-exciton coupled field description. 

The potential V(r) in Eq. ^ describes a defect, which 
can be either naturally present in the cavity mirror [3] 
or it can be created by an additional laser [17] . Later 
on, we will assume the defect to be point-like V(r) = 
gyS(r) and weak, so that we can apply the linear response 
approximation llj. In this treatment, one divides the 
response of the LP field in a mean-field component ipQ 
corresponding to the case when the perturbing potential 
is absent, and a fluctuation part 8ip(r,t) reflecting the 
linear response of the system to the perturbing potential: 



ipLp(r,t) = e 



[e*V r Vo + ^(r,t)] 



(4) 



By substituting Q into ([!]), we obtain a mean-field equa- 
tion and, by retaining only the linear terms in the fluc- 
tuation field and the defect potential, the following first 
order equation in <5?/;(i",i): 



id i 



Sip 



= L 



8ip 
Sip* 



V(r) 



■0 o e ik 



where the operator C is given by: 



r _ ( ^Lp(-iy) - in 



-2ik„ 



(5) 



(6) 




-3 -2-10 1 

Sk x [mc s ] 



2-1012 

Sk x [mc s 



FIG. 1: (Color online) Collective excitation spectra for the 
subsonic (thick solid [black] line at v p — 0.2c s , with c s = 
V 5lV-'o| 2 / m ) and supersonic (dashed [red] line at v p = 1.9c s ) 
regimes and for an interaction-renormalised pump detuning 
A p = -0.3<# | 2 (a, b), A p = (c, d), A p = 0.3 ff |^o| 2 (e, 
f) and A p = 2.3(?|V'o| 2 (g, h). Real parts of the spectra are 
plotted in the left panels and the corresponding imaginary 
parts in the right panels for k — l.lg\ipo\ 2 — note that in our 
description the spectrum imaginary parts do not depend on 
the fluid velocity v p . 



with uj~lp = ^>lp ~ w p + 2.g | ipo \ 2 ■ We are not interested 
here in solving the complex cubic mean-field equation for 
^0) as this has been already widely studied pQ. Rather, 
we want to study the response of the system to the pres- 
ence of the defect and how different behaviours of the 
onset of dissipation can be described in terms of the dif- 
ferent excitation spectra one can get for polaritons reso- 
nantly pumped close to the bottom of the LP dispersion. 



where v p = k p /m is the fluid velocity, and m is the LP 



mass, m 



2m c [l 



A" 



A) 2 



This simplification allows one to describe the complex 
spectrum in terms of three parameters only, namely the 
fluid velocity v p , the interaction-renormalised pump de- 
tuning 



WLP(O) 



2^ 



(9) 



A. Spectrum of collective excitations 



and the LP lifetime k: 



The spectrum of the collective excitations can be ob- 
tained by diagonalising the operator C in the momentum 
space representation: 



C 



-9% 



9^0 
-u) LP (5k - k p 



where, 5k = k — k p . The description of the spectrum 
simplifies in the case when the pumping is close to the 
bottom of the LP dispersion, that can be approximated 
as parabolic 



Lo LP (dk ± k p ) ~ lo lp (Q) 



Sk 2 



2m 2m 



± 5k ■ w, 



(8) 



: (k) = Sk ■ Vp - in ± ^eiSk) [e{Sk) + 2g\ifj 



(10) 



where e(k) = ^— — A p . If energies are measured in units 
of the mean-field energy blue-shift gIV'ol 2 (we will use 
the notation A' p = A p /g\tp \ 2 and k' = K/g\tp \ 2 ), then 
the fluid velocity v p is measured in units of the speed 
of sound c s = \J g\ijjo\ 2 /m. In order to make connection 
with the current experiments, note that, for blue-shifts 
in the range gIV'ol 2 — 0.1 — 1 meV, typical values of the 
speed of sound c s are 0.8 — 2.7 x 10 6 m/s. Similarly, for 
common values of the LP mass, the range in momenta in 
Fig. [l] comes of the order of Sk x ~ 0.2 — 0.8 /im -1 . 
The spectrum ( 10 ) can be classified according to 



the sign of the interaction-renormalised pump detuning 
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A p [IE] — see Fig. [T] For A p < [panels (a,b)], the 
real part of the spectrum is gapped while the imaginary 
part is determined by the polariton lifetime k only. If 
one applies the Landau criterion making reference to the 
real part of the spectrum only, then one finds a critical 
velocity 



1 



I A' 



AilflA' 



2) > 1 



(11) 



always larger than the speed of sound for A p < 0. If 
the fluid velocity is subcritical, v p < v c (see [black] solid 
lines in Fig. [ija)), then no quasiparticles can be excited 
and thus, for infinitely living polaritons k — > 0, the fluid 
would experience no drag when scattering against the 
defect. For supercritical velocities instead, v p > v c see 
[red] dashed lines in Fig.jlja), one expects dissipation in 
the form of radiation of Cherenkov-like waves from the 
defect into the fluid. In the supercritical regime, the set 
of wavevectors k for which 5R[o; + (k)] = form a closed 
curve in the k-space with no singularity of the derivative, 
i.e., in other words, the radiation can be emitted in all 
possible directions around the defect. This, as we will 
see in the next section, will imply that the drag force for 
k — y goes abruptly, rather than continuously, from zero 
at v p < v c to a finite value at v p > v c . 

The spectrum gap closes to zero in the resonant sit- 
uation at A p — 0, when the two branches w ± (k) touch 
at 5k = [panels (c,d) of Fig. [Tj: Here, the real part of 
the spectrum displays the standard linear dispersion at 
small wavevectors as for the weakly interacting bosonic 
gases, with the slope given by c s ± v p . The imaginary 
part, as in the previous case, is constant and equal to 
— k. It is clear therefore that in this case, when n — > 0, 
one recovers the equilibrium results valid for weakly in- 
teracting gases [13 H8] , where the critical velocity for 
superfluidity equals the speed of sound, v c = c s , and the 
drag displays a threshold like behaviour. Here, in the 
supersonic regime v p > v c , the close curve 3?[w + (k)] = 
has instead a singularity, resulting in the standard Mach 
cone of aperture 9, sin# = c s /v p , inside which radiation 
from the defect cannot be emitted [15] . 

Finally, for A p > 0, the real parts of the particle w + (k) 
and hole w~(k) branches of the spectrum touch together 
in either one [A p < 2, see panels (e,f)] or two [A p > 2, see 
panels (g,h)] separate regions in momentum space. In the 
same regions, the corresponding imaginary parts instead 
split. With a somewhat abuse of language, we call these 
kinds of spectrum, diffusive-like. We note that, clearly, 
these spectra have no correspondence in equilibrium sys- 
tems, because a finite polariton lifetime k is needed in or- 
der for these modes to be stable, Q ! [a; ± (k)] < 0. We also 
note that for these spectra, even if considering only the 
real part of the collective excitation spectrum, as soon as 
the fluid is in motion v p > 0, dissipation in the form of 
waves is possible. However, we will see that similarly 
to the case of polaritons non-resonantly pumped [15], 
when decreasing k (and accordingly A p in order to have 
stable solutions), this situation connects continuously to 




FIG. 2: (Color online) Drag force F as a function of the 
fluid velocity v p for different values of the pump detuning 
, = -0.3g|4'o| 2 (a), A p = (b), and A p > (c), and 



for different values of the polariton lifetime 
the notation k — K,/g\ipo\ 2 , A' p = A/g\ipo\ 2 - 



here, we use 



the previous case, where a threshold-like behaviour with 
v c = c s was found. 

We will see in the next section how these different 
spectra imply only two qualitatively different types of 
crossover of the drag force as a function of the fluid ve- 
locity, for either A p < or A p > pump detunings. 



III. DRAG FORCE 

The steady state response of the system to a static and 
weak defect can be evaluated starting from Eq. (j5j): 



y(r)e lk ^ r V'o 
-V{r)e' lk "- r ^ 



For a point-like defect, this can be written in momentum 
space as: 



^ s (k + k p ) 



-gy0o(e(k) - k-Vp + in) 
e(k)[e(k) + 2 3 |^| 2 ]-(k-v p -m)2 ' 



while the other component Sip* (k p — k) can be obtained 
by complex conjugation and by substituting k i— > — k. 
The drag force exerted by the defect on the fluid is given 
by HI]: 



F = -J dv\i, LP {v,t)\ 2 V{V{v)) 



(12) 
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FIG. 3: (Color online) Drag force F as a function of the 
inverse polariton lifetime re' = re/ (g\ ipo\ 2 ) in the (a) subcritical 
regime (v p — 0.2c s ) and (b) supercritical regime (v p — 1.9c s ). 
In both cases we have fixed A p = — 0.3c? ] XE'o 1 2 (v c — 1.46c s ) 
but these results are qualitatively similar for any other value 
of the pump detuning. We plot in the insets the normalised 
real-space wavefunction |'0Lp(r)| 2 /|t/;o| 2 for two specific values 
of re' = 0.2 and re' = 0.8. 



and, in the steady state linear response regime, we ob- 
tain: 



9v 



dk 



ik [ip^Sip s (k + k p ) + ip Q Sip* s (kp - k)] 

dk ike(k) 
(2tt) 2 w+(k)w-(k) 



(13) 



The drag is clearly oriented along the fluid velocity v p , 
i.e., F = Fv p . If k —> 0, then the integral in Eq. ( fl3| is 
finite only if poles exist when 3?[w ± (k)] = 0, i.e., when 
quasiparticles can be excited, in agreement with the Lan- 
dau criterion. For finite polariton lifetimes, however, it 
is clear that the integral will be always different from 
zero for v p > 0. We now analyse the behaviour of the 
drag force as a function of the fluid velocity for the three 
(A p = 0, Ap > 0, and A p < 0) different spectra illus- 
trated in the previous section. 

For the linear spectrum, at A p = 0, in the equilibrium 
limit, k — > 0, we recover for the drag the known result of 
weakly interacting Bose gases in two dimensions (11) : 



with a threshold-like behaviour at a critical fluid veloc- 
ity equal to the speed of sound c s . This limiting re- 
sult is plotted as a bold gray line in the panels (b,c) of 
Fig. [2~| For A p = and finite lifetimes k, we find a 
smooth crossover from the subsonic to the supersonic 
regime, with the drag being closer to the equilibrium 
threshold behaviour for decreasing k (see Fig. [2jb)). A 
finite lifetime tends to increase the value of the drag in 
the subsonic region v p <C v c , giving place to a residual 
drag force, similar to what was found in the numerical 
simulations of Ref. [5 . Instead, in the supersonic region 
tip > ti c , the finite lifetime tends to decrease the value of 



the drag. In the case of diffusive-like spectra at A p > 
the situation is qualitatively very similar to the resonant 
case (see Fig. [2jc)) , with the difference that now, in or- 
der to have stable solutions, we can decrease the value 
of the lifetime only by decreasing accordingly also the 
value of the pump detuning A p . The crossover for both 
Ap = and A p > is also qualitatively very similar to 
the case of non-resonantly pumped polaritons [TS] , where 
the spectrum of excitation is in that case diffusive-like. 

In the case of gapped spectra, the situation is however 
qualitatively different (see Fig. [2^ a)). For infinitely living 
polaritons, k — > 0, the drag force can also be evaluated 
analytically and its expression is similar to Eq. ( 14 ) , but 



with a critical velocity larger than the speed of sound, 
which expression is given in Eq. (Ill: 



F 

(mc s ) 3 gv/g 



v p /c s 



1 



e(v p - v c ) 



(15) 



Therefore now the drag experiences a jump for v p = v c , 
rather than a continuous threshold as for the resonant 
case Ap = 0. As already mentioned in the previous 
section, this discontinuous behaviour of the drag for the 
gapped spectra is connected to the fact that, as soon as 
quasiparticles can be excited by the defect at v p > v c , 
Cherenkov-like waves can be immediately emitted in all 
directions, rather than being restricted in a region out- 
side the Mach cone like before. For A p = 0, the cone was 
gradually closing with increasing the fluid velocity. 

Both the increase of the value of the drag in the sub- 
critical region as a function of the polariton lifetime and 
the decrease in the supercritical region, are behaviours 
common to all the types of spectra. We plot the drag 
force as a function of re in Fig. [3j for two values of the 
fluid velocity v p and a specific value of the pump detun- 
ing Ap, though we have checked that the following results 
are generic. For v p < v c , we find that the residual drag 
is a finite-lifetime effect only, and, in agreement with the 
results of Ref. [5], we find that, well below the critical 
velocity, the drag force goes linearly to zero for k 0. 
In the resonant case A p = 0, the slope of the drag for 
v p <C c s can be evaluated analytically starting from the 



expression (131 



F 



(v p /c s ) 2 



v p /c s 



-0(v P 



(14) 
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2c, 



(mc s ) 3 g%r/g k^o ttv p \ ^1 - {v p /c s )'- 
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The residual drag in the subsonic regime is an effect of the 
broadening of the quasi-particles energies: Even when the 
spectrum real part does not allow any scattering against 
the defect (e.g., for A p < 0), the broadening produces 
some scattering close to the defect. This results in a per- 
turbation of the fluid around the defect, asymmetric in 
the direction of the fluid velocity (see panel (a) of Fig. [3]), 
similarly to what it was obtained in Ref. [5! . Instead, in 
the supersonic regime, the drag force is weaker in the non- 
equilibrium case respect to the equilibrium one. This is 
caused by the finite lifetime tending to suppress the prop- 
agation of the Cherenkov waves away from the defect, as 
shown in panel (b) of Fig. [3j 



IV. CONCLUSIONS AND DISCUSSION 

To conclude, we have analysed the linear response to 
a weak defect of resonantly pumped polaritons in the 
pump-only state and we have been able to determine 
two different kinds of threshold like behaviours for the 
drag force as a function of the fluid velocity. In the 
case of either zero or positive pump detuning, one can 
continuously connect to the case of equilibrium weakly 
interacting gases, where the drag displays a continuous 
threshold with a critical velocity equal to the speed of 
sound. However, for negative pump detuning, where the 
spectrum of excitations is gapped, the drag shows a dis- 
continuity with a critical velocity larger than the speed of 
sound. In this sense, the case of coherently driven micro- 
cavity polaritons in the pump-only configuration displays 
a richer phenomenology than the case of polariton super- 
fluids non resonantly pumped. It would be interesting to 
perform a similar analysis in the case of polaritons in the 
optical parametric oscillator regime, where polaritons are 
parametrically scattered from the pump state to the sig- 
nal and idler states. Here, the spectrum of excitations 
has been already determined in Ref. [19j , however it is 
far from clear what are the conditions for subcritical, 
superfluid, behaviour in a fluid characterised by three 
distinct currents, and how the link between signal and 
idler imposed by the parametric scattering influences the 
scattering of both fluids against a defect. 
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Appendix A: Gross-Pitaevskii equation for the lower 
polariton field 

If one starts from a descriptions of polaritons in terms 
of separate exciton and cavity photon fields, a rota- 
tion into the lower and upper polariton basis, followed 
by neglecting the occupancy of the upper polariton 
branch, results in the following Gross-Pitaevskii equation 
for the lower polariton (LP) field in momentum space 
^p(r,t) = E k e ikr ^P,k(t) [18: 

id t ipLPM = f P e^ pt S^ kp + [ui L p{k) - in(k)} ?/>LP,k+ 

5 , k,k 1 ,k 2 V'LP,k 1 +k 2 -k ? /'LP,k 1 V'LP,k 2 + 



ki,k 2 



Sfe VVk-kiV'LP^Sfci , (Al) 



k 



where fc(fe) = nxc\ + k-cs\ is the effective LP decay 
rate, 5 k ,ki,k 2 = gxCfcC| kl+ k 2 -k|Cfe 1 c/ C2 is the interaction 
strength, and where V(r) = ^ k e lkr Vk- In these expres- 
sions, the coefficients 



2 2 

c ki s k 



1± 



X 



(A2) 



are the Hopfield coefficients used to diagonalise the free 
polariton Hamiltonian. We want here to justify the sim- 
plified description done in Eq. ([I]). If we follow the linear 
response expansion as in (El), the operator £ in momen- 
tum space analogous to (I7J) reads as: 



^ fu} LP (5k + k p ) -in(5k + k p ) ffxCfe p C 5k+kp c 5k _ kj ,Vo \ (A3) 
k,kp 1 -gxcLQ k+ k p c<5k- kp V'c> 2 -ULp(5k-k p ) - in(Sk-kp) J ' 



where now ujLp(5k ± k p ) = ujLp(Sk ± k p ) 



-5k±k, 



IV'ol • It is easy to show that the eigenval- of 



ues of this operator coincide with our approximated ex- 
pressions ( 10 1 in the limit of 5k <C k p , when c^ k±k ~ c\ , 



, (5k±k„ 



sj, and when we can simply rename g = gx c 



and k = n{kp). It is interesting to note that, even if we 



would retain the linear terms in k p ■ 5k in the expansion 



'<5k±k„ 



this would result in a rcnormalisation of the 



fluid velocity v p in the expression ( 10 1 which takes into 
account the blue-shift of the lower polariton dispersion 
due to the interaction. 
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